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HE. 【目的 ] MiRNAs 在 昆虫 变态 发 育 过 程 中 发 挥 非常 重要 的 作用 。 对 家 彼 Bombyx mori miRNAs 
pi O miRNAs A457 RRRSRAVD TM, [7A] ERE 5 BA 
2 KOS MoH CIES RRMA 20E 后 ,qRT-PCR 检测 miR-2769 在 家 看 脂肪 体 中 的 表达 ;通过 生物 
信息 学 方法 预测 家 和 看 miR-2769 的 鞠 基 因 ; 利 用 双 荧 光 酶 报告 载体 系统 分 析 miR-2769 与 预测 靶 基 
因 BmE75B 的 互 作 ;qRT-PCR 检测 miR-2769 及 其 靶 基 因 BmE75 不 同 剪接 体 在 家 蚕 不 同 发 育 时 期 
(幼虫 \ 肾 和 成 虫 ) 和 幼虫 不 同 组 织 ( 头 、 表 皮 、 丝 腺 、 脂 肪 体 、 精 烛 、 卵 集 、 蕊 民 管 .中 肠 和 血 淋 巴 ) 中 
的 表达 量 。【 结 果 】 研 究 结果 表明 ,miR-2769 可 通过 与 家 看 BmE75B 的 3'UTR 区 结合 位 点 的 互 作 ， 
显著 抑制 荧光 素 酶 报告 基因 的 表达 。qRT-PCR 结果 表明 ,miR-2769 和 BmE75A/BmE75B 在 20E 诱 
导 家 蚕 脂肪 体 中 表达 趋势 相反 。 时 空 表达 分 析 结 果 表 明 ,miR-2769 与 BmE75 HARB BEARER 
甜 不 同 发 育 时 期 和 不 同 组 织 中 均 具有 特异 性 表达 特征 。 在 家 和 剧变 态 发 育 的 不 同 阶段 ,miR-2769 和 
BmE754 的 表达 量 均 较 低 ,而 BmE75B 在 肾 期 表达 量 较 高 BmE75C 在 5 龄 未 幼虫 和 肾 早 期 有 极 高 
表达 水 平 。 此 外 ,miR-2769 与 BmE75 不 同 剪 接 体 均 存在 一 定 程度 的 表达 负 相 关 。 在 家 看 5 Hoh R 
血 淋 巴 中 miR-2769 可 促进 BmE754 的 表达 ,在 脂肪 体 中 miR-2769 可 促进 BmE75B 的 表达 ,而 在 其 
他 组 织 中 miR-2769 抑制 BMETS 不 同 剪 接 体 的 表达 。 【结论 】] Ke miR-2769 可 通过 与 BmE7SB 的 
3'UTR 区 i 互 作 对 BmETS 不 同 剪接 体 的 表达 进行 负 ae 
关键 词 : KR; 变态 发 育 ; SRR; miRNAs; 靶 基 因 ; qRT-PCR; 表达 分 析 
中 图 分 类 号 : Q966 文献 标识 码 : A 。 ”文章 编号 : 0454-6296 (2020)03-0245-10 
Identification and expression analysis of miR-2769 target gene in Bombyx 
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Abstract; [Aim] MiRNAs play very important roles in insect metamorphosis. The study of miRNAs and 
their target genes in Bombyx mori will help to elucidate the molecular mechanism of miRNAs involved in 
the metamorphosis of this insect. [Methods] The expression of miR-2769 in the fat body of B. mori was 
detected by qRT-PCR after injecting molting hormone, 20-hydroxyecdysone (20E) , into the hemolymph 
of the day-2 5th instar larvae. The target genes of miR-2769 of B. mori were predicted by bioinformatics 
methods. The interaction between miR-2769 and its predicted target gene BmE75B was analyzed by dual 
luciferase reporter assay system. The expression levels of miR-2769 and different splice isoforms of its 


target gene BmE75 at different developmental stages (larva, pupa and adult) and larval tissues (head, 
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epidermis, silk gland, fat body, testis, ovary, Malpighian tubules, midgut and hemolymph) of B. mori 
were detected by qRT-PCR. [Results] The results showed that miR-2769 could significantly inhibit the 
expression of luciferase reporter gene by interacting with the 3’UTR binding site of BmE75B of B. mori. 
The qRT-PCR results showed that miR-2769 and BmE75A/BmE75B presented the opposite expression 
trends in the fat body of B. mori post induction by 20E. The spatio-temporal expression analysis showed 
that miR-2769 and different splice isoforms of the target gene BmE75 had specific expression 
characteristics in different developmental stages and larval tissues of B. mori. The expression levels of 
miR-2769 and BmE75A were low in B. mori at different developmental stages; however, the expression 
level of BmE75B was relatively high in the pupal stage and that of BmE75C was extremely high in the late 
5th instar larvae and early pupae. Furthermore, the expression of miR-2769 was negatively correlated 
with that of different splice isoforms of BmE75. miR-2769 promoted the expression of BmE75A in the 
hemolymph and BmE75B in the fat body of the Sth instar larvae of B. mori, while it inhibited the 
expression of different splice isoforms of BmE75 in other tissues. [ Conclusion] MiR-2769 can negatively 
regulate the expression of different splice isoforms of BmE75 by interacting with the 3’ UTR region of 
E75B in B. mori. 

Key words: Bombyx mori; metamorphosis; molting hormone; miRNAs; target gene; qRT-PCR; 


expression analysis 





MicroRNAs ( miRNAs ) 是 一 种 内 源 性 小 分 子 非 
编码 RNA ,长 度 约 22 ni, 广泛 存在 于 原生 动物 、 病 
ae HH Al oh wy R, 且 物 种 间 具 有 高 度 保 守 人 性 
(Ambros, 2004; Bartel, 2004), miRNAs 通过 与 靶 
基因 的 非 翻 译 区 或 编码 区 结合 ,在 转录 后 或 翻译 水 
平 上 调节 靶 基 因 mRNA 的 降解 或 蛋白 的 表达 ,参与 
调控 细胞 增殖 分 化 (Li et al., 2013) YT (He et al., 
2013; Li et al., 2018). A We ( Nelson et al., 2014; 
Texada et al., 2019; Zhao et al., 2019) , 4% (Silwal 
et al., 2019; Xu et al., 2019) 、 发 育 和 代谢 (He and 
Hannon, 2004; Ameres and Zamore, 2013 ) 等 几乎 所 
有 的 生命 活动 过 程 。lin4 和 lin-7 是 最 早 在 秀丽 隐 
杆 线 虫 Caenorhabditis elegans 中 发 现 的 miRNAs( Lee 
et al., 1993) ;在 哺乳 动物 中 , 约 50% 的 编码 蛋白 基 
因 表 达 受 miRNAs 的 调控 ( Krol et al., 2010); 
miRNAs 同样 在 昆虫 中 广泛 存在 ,并 在 昆虫 变态 发 
育 和 生殖 过 程 中 发 挥 重 要 作用 (Behura，2007 ; 
Lucas and Raikhel, 2013; Liu et al., 2018, 2019) 。 

BE HAAS AS Ac A EH WE BZ HA ZR 20E 和 保 幼 激素 





























,细胞 核 激 素 受 体 E75 (nuclear hormone receptor 
E75) 基因 最 早 在 黑 腹 果 晶 Drosophila melanogaster 
中 被 鉴定 ,该 基因 在 果 蝇 中 通过 3 ' 端 外 显 子 的 可 变 
剪接 编码 DmE75A, DmE75B, DmE75C 和 DmE75D 
4 个 异 构 体 ,其 中 DmE75A 和 DmE75C 的 DNA 结合 
结构 域 含 有 2 个 锌 指 结构 ,DmE75B 含有 1 个 锌 指 
结构 ,DmE75D 没有 DNA 结合 结构 域 (Feigl et al., 
1989; Segraves and Hogness, 1990; Bernardo et al., 
2009), EA HEP ETS 基因 具有 BmE75A, BmE75B 
和 BmE75C 3 个 可 变 剪接 体 , 且 其 编码 蛋白 同样 具 
有 相同 的 C 末端 和 不 同 的 N 末端 (Segraves and 
Hogness, 1990), 对 BmE75A, BmE75B 和 BmE75C 
EZKE BmN 细胞 .前 胸腺 和 脂肪 体 中 的 表达 研究 结 
果 也 表明 ,BmE75A4，BmE75B 和 BmE75C 具有 不 同 
的 组 织 和 时 期 特异 性 , 而且 三 者 的 表达 与 20E 的 滴 
度 密切 相关 ( Li et al., 2015) 。 近 年 研究 表明 ,E75 FE 
因 可 参与 不 同 昆 虫 幼虫 到 贤 的 变态 发 育 过 程 (Guo 
et al., 2016; Liu et al., 2018; Peng et al., 2019) 。 
Zi Bombyx mori 作为 重要 的 经 济 昆虫 ,是 研 

























































































(juvenile hormone, JH) HEJJ IE ,研究 发 现 ,20E 不 
仅 可 抑制 PIBK-Akt 信和 号 通路 而 促进 自 噬 的 起 始 , 同 
时 可 与 20E 受 体 EcR-USP 的 异 源 二 聚 体 结合 ,启动 
20E 初级 应 答 基 因 Br-C, E74, E75 和 E93 的 表达 ， 
进而 通过 与 特异 DNA 序列 结合 启动 晚期 应 答 基 因 
的 表达 (Riddiford et al., 2000; Yin and Thummel, 
2005)。 作 为 20E 信号 通路 中 的 一 个 初级 应 答 基 


















































究 昆 虫 变态 发 育 分 子 机 制 的 重要 模式 生物 。 家 看 的 
变态 发 育 过 程 由 20E 和 JH 协同 调控 。 在 昆虫 变态 
发 育 过 程 中 需 经 历 复杂 的 器 官 重建 过 程 ,脂肪 体 的 
日 哈 和 凋 亡 对 家 看 变态 发 育 过 程 十 分 重要 ( Ihry and 
Bashirullah, 2014; Hu et al., 2016)。 本 研究 通过 生 
物 信 息 学 方法 分 析 预 测 家 看 miR-2769 在 20E 处 理 
后 家 不 脂肪 体 中 差异 表达 基因 中 的 关联 可 基因 , 利 
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MIKIR BET Fa SE AL fr Ul AA BE XT miR-2769 与 其 
WEA BmE75B 的 互 作 进 行 鉴定 ,进一步 通过 qRT- 
PCR 检测 miR-2769 、 靶 基因 BmE75B 及 其 不 同 剪接 
体 BmE75A 和 BmE75C 在 20E 处 理 后 家 春 脂 肪 体 及 
不 同 发 育 时 期 和 不 同 组 织 的 表达 特征 , 从 而 阐明 
miR-2769 及 其 靶 基 因 BmE75B Ee RRA Bt 
程 中 的 作用 。 




















1 材料 与 方法 


1.1 材料 与 试剂 

供 试 家 春 为 二 化 性 的 大 造 p50 品系 ,常规 条 件 
下 对 乔 儿 进行 俊 青 ,幼虫 用 新 鲜 桑 叶 饲 养 。 收 集 5 
HER 2 天 家 乔 幼 虫 进 行 血 淋 巴 注射 晓 皮 激素 20E 
(2 pe/wh) ,每 头 家 和 蛋 注 射 5 ng, 对照 注 射 等 体积 
菌 水 稀释 10 倍 的 无 水 乙醇 溶剂 。 注 射 后 2, 6, 12 
和 24 h 后 分 别 取 处 理 和 对 照 的 家 蛋 幼 虫 脂肪 体 。 
家 看 头 、 表 皮 、 丝 腺 、 脂 肪 体 、 精 人 巢 、 卵 集 、 马 氏 管 中 
肠 和 血 淋 巴 的 组 织 取 自 正常 家 在 5 龄 第 3 天 幼虫 ; 
在 4 龄 幼虫 .5 龄 幼虫 、 肾 期 至 羽化 后 3 d 成 虫 期 间 
每 天 分 别 随 机 取 整 虫 ,幼虫 去 中 肠 。 每 个 样品 3 H, 
每 组 来 自 5 头 个 体 。 所 有 样品 液 氮 速冻 后 - 80°C vk 
箱 保 存 备用 。 

大 上 肠 杆菌 Escherichia coli Top10 细胞 ,人 上 胚胎 肾 
细胞 HEK293T 细胞 ,psiCheck2 载体 为 实验 室 保存 ， 
用 于 基因 克隆 的 载体 和 试剂 有 pMD-19T 载体 
(TaKaRa, Ki) ,Trizol Reagent( Invitrogen, 美国 )， 
RevertAid First Strand cDNA Synthesis Kit 
(Thermofisher ,美国 ) ,限制 性 内 切 酶 (Fermentas , 美 
E), 质粒 DNA 提取 试剂 盒 和 胶 回 收 试剂 盒 
(Axygen, 美国 ) , Taq 酶 (Vazyme， 南 京 ) ,T4 DNA 
连接 酶 (Thermofisher, 美国 ) ;用 于 细胞 培养 的 试剂 
有 DMEM 高 糖 培养 基 ( Hyclone, 美国 ) , 胎 牛 血清 
(BI， 以 色 列 ), 细胞 转 染 试剂 为 X-tremeGENE 
siRNA Transfection Reagent (Roche, mE), M356 
素 酶 报告 基因 检测 试剂 盒 ( Promega, 美国 ) , qRT- 
PCR 试剂 为 FastStart 通用 型 SYBR® Green 预 混 液 
(Roche, 瑞士 ) 和 All-in-one™ miRNA qRT-PCR 检 
测试 剂 盒 (Cenecopia, 美国 ) 。 

1.2 R24 RNA 的 提取 

将 1.1 PARARE lh OP Sl) RTS, PY 
100 mg 粉末 加 入 1 mL Trizol 试剂 ,按照 RNA 提取 
说 明 书 进行 抽 提 ,RNA 沉淀 充分 溶解 后 用 超 微量 核 
酸 蛋 白 测 定 仪 (NanoDrop ND-2000) 对 总 RNA 的 浓 




































































度 和 纯度 进行 测定 。 
1.3 RE miR-2769 的 序列 和 表达 分 析 

根据 前 期 测序 获得 的 成 熟 miR-2769 的 序列 
( AUAUAUUAUCAGAUUUUCGGUC ) 设计 上 游 引 物 
(5'-ATATATTATCAGATTITTCGGTC-3’) ,以 All-in-one™ 
miRNA qRT-PCR 检测 试剂 盒 中 的 通用 的 接头 引物 
(Genecopia, 美国 ) 为 下 游 引 物 ,对 1.2 节 家 乔 脂 肪 
体 RNA 进行 反 转 录 , 以 反 转 录 合 成 的 cDNA 为 模 
板 , 利 用 All-in-one™ miRNA qRT-PCR 检测 试剂 盒 
进行 PCR 扩 增 , 扩 增 条 件 : 预 变 性 95%C 5 min; 
95° 15 s, 55%C 30 s, 72% 30 s, 共 30 个 循环 ;最 后 
72% 延伸 10 min。 扩 增产 物 纯化 后 连接 pMD-19T 
载体 ,连接 产物 转化 大 肠 杆菌 Top10 感受 态 细胞 , 挑 
取 鉴 定 后 的 阳性 克隆 送 深 圳 华 大 基因 科技 有 限 公 司 
进行 测序 。 

利用 All-in-one™ miRNA qRT-PCR 检测 试剂 盒 
对 miR-2769 在 1.1 节 中 20E 处 理 组 与 对 照 组 的 2 ， 
6, 12 和 24 h 的 脂肪 体 组织 中 的 表达 进行 qRT-PCR 
检测 。 内 参 基因 为 家 看 5S rRNA (GenBank 登录 号 : 
M35394. 1), 引物 序列 为 SS-F: 5'-CGTCCGATCA 
CCGAAGTCAAG-3'; 5S-R: 5'-AGGCGGTCACCCAT 
CCAAGTA-3'。 扩 增 条 件 : 95°C 5 min; 95%C 15 s, 
55°C 30 s, 72C 30 s, 共 40 次 循环 ;最 后 72% 延伸 
10 min ,溶解 曲线 范围 为 65Y ~9S% 。 

1.4 RE miR-2769 靶 基 因 的 预测 及 表达 分 析 

课题 组 前 期 对 20E 处 理 前 后 的 家 不 脂肪 体 进 
行 了 转录 组 测序 , 据 此 利用 软件 miRanda 和 生物 信 
息 学 方法 对 20E 处 理 前 后 家 和 蛋 脂 肪 体 转录 组 数据 
中 差异 表达 基因 中 的 miR-2769 4E SEA ETT THD 
利用 ClustalW2 在 线 分 析 软 件 (https: // www. ebi. ac. 
uk/Tools/msa/clustalw2/) 对 miR-2769 和 靶 基 
BmE75B( GenBank 登录 号 : NM_001112610. 1) 及 其 
不 同 剪接 体 BmE75A ( GenBank 登录 号 : NML 
001112609.1) 和 BmE75C ( GenBank 登录 号 : NM_ 
001043577. 1) 的 序列 进行 比 对 分 析 。 

以 1.2 节 提取 的 家 看 幼虫 脂肪 体 RNA 为 模板 ， 
利用 RevertAid First Strand cDNA Synthesis Kit 反 转 
录 试 剂 盒 合成 cDNA, Fl] FA FastStart 通用 型 SYBR® 
Green 预 混 液 对 miR-2769 靶 基 因 BmE75B 及 其 不 同 
剪接 体 BmE75A 和 BmE75C 进行 qRT-PCR 检测 。 
扩 增 条 件 : 95°C 5 min; 95° 15 s, 55% 30 s, 72% 
30 s, 共 40 次 循环 ;最 后 72% 延伸 10 min ,溶解 曲线 
范围 为 65%C ~95°C., BmE75B 的 引物 序列 为 BmE75B- 
F: 5'-CACTCGGTGCTAGTGAGCAT-3'; BmE75B-R; 
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5'-TGCAGGGCCGGTATTGTATC-3’, BmE75A 的 引 
物 序列 为 BmE75A-F: 5'-GGAACGTGAACCTGAA 
TTGC-3'; BmE75A-R; 5'-TGCAGGGCCGGTATTGT 
ATC-3', BmE75C 的 引物 序列 为 BmE7SC-F; 5'- 
AGTGTTATCCGAAGTTGTCAC-3'; BmE75C-R: 5'- 
GAATGATGGAGGCCTTGAGCA-3'。 上 述 引 物 序 列 分 
别 对 应 于 BmE75B, BmE75A 和 BmE75C 的 N 端 特 
RRS, NBA NARA Bmactin43( GenBank 登录 
号 : U49854 ) ,引物 序列 BmA3-F; 5'-ATGTGCGA 
CGAAGAAGTTGC-3'; BmA3-R; 5’-GTCTCCTACGTA 
CGAGTCCT-3’, 
1.5 RÆ miR-2769 的 合成 和 BmE75B 3'UTR X 
荧光 素 酶 报告 载体 构建 

根据 成 熟 miR-2769 的 序列 及 其 与 预测 靶 基 因 
BmE75B 的 互 作 结合 位 点 序列 ,由 上 海 吉 玛 生 物 技 
术 有 限 公司 设 计 并 合成 miR-2769 的 野生 型 模拟 物 
miR-2769 mimics (5’-AUAUAUUAUCAGAUUUUCGG 
UC-3') 互 作 位 点 突变 后 的 miR-2769 mutant (5'- 
AAUAUAAUUCAGAUUUUCGGUC-3') 及 阴性 对 照 
miR-2769 NC(5'-UUCUCCGAACGUGUCACGUTT-3’) 。 

LX RevertAid First Strand cDNA Synthesis Kit 反 
转录 试剂 盒 合成 的 家 看 脂肪 体 cDNA 为 模板 ,利用 
BmE75B 3'UTR 序列 的 特异 引物 ,PCR 扩 增 277 bp 
目的 片段 ,纯化 后 扩 增 产物 和 psiCheck2 载体 分 别 
用 限制 性 内 切 酶 Xho IAI Not I 进行 双 酶 切 , 酶 切 产 
物 切 胶 纯 化 后 用 T4 DNA 连接 酶 连接 ,连接 产物 转 
化 大 肠 杆菌 Top10 感受 态 细胞 ,经 菌 液 PCR 和 双 酶 
切 鉴定 后 的 psiCheck2/BmE75B 3'UTR 阳性 克隆 送 
深圳 华 大 基因 科技 有 限 公 司 进 行 测序 验证 。 引 物 序 
列 为 BmE75B 3' UTR-Xho | -F; 5'-AGCTCGAGGC 
CCCGACCTACTITAAA-3'; BmE75B 3'-UTR-Not | - 
R: 5'-ATGCGGCCGCAATTICGACAACCTGTC-3’ ( F 
划 线 为 酶 切 位 点 )。 
1.6 RE miR-2769 与 靶 基 因 BmE75B 的 互 作 分 析 

HEK293T 细胞 于 37C ,5% CO, 的 条 件 下 ,用 
10% 胎 牛 血清 的 DMEM 高 糖 培养 基 培 养 。 将 生长 
状态 良好 的 细胞 提前 d 铺 于 24 孔 板 中 用 于 转 染 
实验 。 将 合成 的 miR-2769 mimics, miR-2769 mutant 
和 miR-2769 NC 分 别 与 psiCheck2/BmE75B 3' UTR 
共 转 染 HEK293T 细胞 , psiCheck2/BmE75B 3'UTR 
质粒 单 转 为 对 照 ,根据 X-tremeGENE siRNA 转 染 试 
剂 的 说 明 书 ,优化 后 每 个 孔 中 质粒 DNA 的 量 为 100 
ng,miRNA 的 终 浓度 为 0.1 pmol/L, 每 组 转 染 3 个 
重复 ,实验 重复 3 次 。 转 染 后 48 h 收集 细胞 ,根据 


























双 荧 光 素 酶 报告 基因 检测 试剂 盒 说 明 书 ,对 细胞 进 
行 裂解 ,利用 化 学 发 光 仪 依次 检测 海 肾 荧光 素 酶 活 
Pe Se ke HEH MT VE, Z W HEK Fee KH 
光 素 酶 的 相对 活性 ( 杨 娟 娟 等 , 2018 ) 。 
1.7 ZE miR-2769 与 靶 基 因 BmE75 不 同 剪 接 体 
的 表达 分 析 

以 1.2 节 提 取 的 家 看 5 龄 幼虫 不 同 组 织 、4 龄 
幼虫 到 3 日 龄 成 虫 整 虫 的 RNA 为 模板 ,利用 All-in- 
one™ miRNA qRT-PCR 检测 试剂 盒 进行 反 转 录 后 对 
miR-2769 进行 qRT-PCR 检测 ,内 参 基因 为 家 看 5S 
rRNA( 方 法 同 1.3 节 )。 同 时 ,利用 RevertAid First 
Strand cDNA Synthesis Kit 反 转 录 试 剂 盒 合成 上 述 
RNA 的 cDNA, ,利用 FastStart 通用 型 SYBR® Green 
预 混 液 对 BmE75B 及 其 不 同 剪接 体 BmE75A 和 
BmE75C 进行 qRT-PCR 检测 , ABE AW RH 
BmactinA3 (方法 同 1.4 节 )。 
1.8 数据 分 析 

双 奖 光 素 酶 活性 检测 结果 ,根据 董 火 虫 荧光 素 
酶 活性 进行 校正 ,以 海 肾 荧光 素 酶 活性 与 莉 火 虫 菊 
光 素 酶 活性 的 比值 为 萤火虫 奕 光 素 酶 的 相对 活性 ， 
利用 Prism 5. 0 作 图 软件 进行 作 图 和 统计 分 析 。 
qRT-PCR 扩 增 结果 用 2-^%" 法 进行 数据 分 析 , 针 对 
不 同 内 参 基因 的 qRT-PCR 扩 增 结果 , 先 通过 计算 两 
个 内 参 基因 Ct 值 的 几何 平均 数 对 不 同 内 参 基因 的 
扩 增 进行 校正 ,再 用 2-^^" 法 进行 数据 分 析 , 最 后 
用 Prism 5.0 软件 进行 作 图 和 统计 分 析 。 


2 结果 


2.1 RE miR-2769 的 序列 和 表达 分 析 

以 家 看 脂肪 体 组 织 RNA 反 转 录 后 的 cDNA 为 
模板 ,PCR 扩 增 miR-2769 , 扩 增 产物 连接 pMD-19T 
克隆 载体 ,连接 产物 经 初步 鉴定 后 进行 测序 ,测序 结 
果 表 明 扩 增 后 miRNA 序列 与 成 熟 miR-2769 序列 完 
全 一 致 。 利 用 qRT-PCR 检测 miR-2769 在 20E 处 理 
后 2, 6, 12 和 24 h 的 家 看 5 龄 第 2 天 幼虫 脂肪 体 
组 织 中 的 表达 情况 。 结 果 如 图 1 所 示 ,与 对 照相 比 ， 
miR-2769 在 20E 处 理 后 2, 6 和 12 h 表达 显著 下 调 
(P<0.05),24 h 时 差异 不 显著 (P>0.05)。 
2.2 RB miR-2769 靶 基 因 的 预测 及 表达 分 析 
通过 生物 信息 学 预测 ,发 现 miR-2769 的 一 个 靶 
基因 是 家 春 的 核 激 素 受 体 基因 BmETS 的 一 个 可 变 
BY BAK BmE75B, miR-2769 第 2 -8 位 种 子 序列 与 
预测 靶 基 因 BmE75B 3' UTR 存在 互补 配对 (图 2: 
A)。 利 用 ClustalW2 在 线 分 析 软 件 对 BmE75B 不 同 
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图 1 qRT-PCR 检测 miR-2769 在 20E 处 理 后 不 同 
时 间 家 看 幼虫 脂肪 体 中 的 表达 模式 


Expression profiles of miR-2769 in larval fat 














Fig. 1 
body of Bombyx mori at different time after 
treatment with 20E by qRT-PCR 
处 理 组 中 往 5 龄 第 2 天 家 和 蛋 幼 虫 血 淋巴 中 注射 赔 皮 激素 20E 
(2 kg/kL) ,每 头 注射 5 pg, 对 照 组 注射 等 体积 无 菌 水 稀释 10 售 
的 无 水 乙醇 溶剂 。 图 3 同 。In the treatment group, 2 pg/pL o 
molting hormone 20E (5 yg per individual) was injected into the 


























hemolymph of the day-2 5th instar larvae, while in the control group 
(CK), an equal volume of ethanol in 10-fold dilution with distille 
water was injected. The same for Fig. 3. CK; 无 水 乙醇 Ethanol. 图 
中 数据 为 平均 值 上 标准 误 , 用 双 尾 配对 上 检验 进行 统计 分 析 。 
Data in the figure are mean + SE, and statistics analysis was performe 
by two tailed, paired t-test. * P <0.05; “无 显著 差异 No significant 
difference (P >0.05). 






































A BmE75B 3'UTR 


剪接 体 的 氨基 酸 序列 比 对 , 分析 结果 表明 ， 
BmE75A, BmE75B 和 BmE75C 的 氨基 酸 序列 具有 
不 同 的 N 端 序列 和 完全 相同 的 C 末端 (图 2: B)。 
同时 BmE75A, BmE75B 和 BmE75C 的 3'UTR 区 
DNA 序列 也 完全 一 致 ,因此 BmE75A 和 BmE75C 的 
3'UTR 区 也 存在 miR-2769 的 结合 位 点 。 

为 了 比较 预测 靶 基 因 BmE7SB 及 其 剪接 体 与 
miR-2769 在 20E 处 理 后 家 不 幼虫 脂肪 体 组 织 中 的 
表达 ,利用 qRT-PCR 检测 BmE75B 在 20E 处 理 组 与 
对 照 组 在 2, 6, 12 和 24 h 脂肪 体 组 织 中 的 表达 。 
结果 如 图 3(A) 所 示 , 与 miR-2769 的 表达 (图 1) 相 
反 ,BmE75B 在 20E 处 理 后 2, 6 和 12 h 表达 显著 上 
调 (P <0. 001) ,24 h 时 表达 差异 不 显著 (P > 
0.05)。 同 时 也 检测 了 BmE75A 和 BmE75C 在 相同 
样品 中 的 表达 情况 ,结果 如 图 3(B) 和 3(C) 所 示 , 与 
BmE75B 类 似 ,BmE754 在 20E 处 理 后 2,6 和 12h 
家 和 看 幼虫 脂肪 体 中 表达 也 显著 升 高 (P <0.001)， 
24 h 时 表达 差异 不 明显 (P > 0.05 ) ; BmE75C 仅 在 
20E 处 理 后 6 h 时 显著 上 调 (P <0.01),2, 12 和 24 
h 时 与 对 照相 比 无 显著 差异 (P >0.05)。 由 此 可 
见 ,在 脂肪 体 中 miR-2769 与 BmE75A 和 BmE75B 的 
表达 存在 负 相 关 , 推 测 在 家 泰 脂肪 体 中 miR-2769 可 
通过 与 靶 基 因 BmE75 3'UTR 的 互 作 调 控 BmE754 和 
BmE75B 的 表达 。 




























































































yo ol G UUU Ay 
GCOGG GGU UG GUGUGT 
= UCU UUA AC  UAUAUA 
yc G UAT T 
miR-2769 
B 
并 MVRTMSCGAELRE----—-—— 
BMETSA: =- SSeS SSS SSS SSS 一 一 -一 MSPDSSYGRYD 
BmE75C MQCYPKLSPKREPPEGLYE I EMLPGARRLELPAPPGKEFRAPVLLAGPSLAPTHSVIQCM 
BmET59B- =a RHSVLVSML 一 一 一 一 一 BARRESSDSGCSSDDSSDVERDSTKESCDR 
BmE75A VSS PTSVDHSLMSSMHKEREPELHIEFDGTTVLCRVCGDKASGFHYGVHSCEGCK 
BmE75C RPPPPPPPPPPPRLLKPPSFEEPSSS IPDLEFDGTTVLCRVCGDKASGFHY GVHSCEGCK 
BmE75B GFFRRS LQQK TQYRPCTKNQQCS TLR INRNRCQYCRLKKCTAVGMSRDAVRFGRVPKREK 
BmE75A GFFRRS LQQK1QYRPCTKNQQCS ILRINRNRCQYCRLKKCIAVGMSRDAVRFGRVPKREK 
BmE75C GFFRRS [QQK [QYRPCTKNQQCS I LR INRNRCQYCRLKKCIAVGMSRDAVRFGRVPKREK 


图 2 miR-2769 FUSED FHM (A) 和 靶 基 因 及 其 不 同 剪接 体 的 序列 比 对 分 析 (B) 


Fig. 2 Prediction of miR-2769 target gene (A) and multiple sequence alignment of the target gene 





and its different splice isoforms (B) 


A: miR-2769 在 靶 基 因 BmE75B 3'UTR 上 的 结合 位 点 预测 Predicted binding site of miR-2769 in 3’UTR of the target gene BmE75B,; B; 靶 基 因 不 同 
剪接 体 的 N 端 氨基 酸 序列 分 析 N-terminal amino acid sequence analysis of different splice isoforms of the target gene. 














250 


昆虫 学 报 Acta Entomologica Sinica 


63 卷 





BmE75BRAXt KI 
BmE75B relative expression level 





2 6 12 24 


处 理 时 间 (h) 


Treatment time 


a 
ied N 
己 in 


一 
wn 


0.5 


BmE75C 相 对 表 
BmE75C relative expression level 





12 24 
处 理 时 间 () 


Treatment time 


2 6 


w 


300 


200 


100 


BE754 相 对 表达 量 


BmE75A relative expression level 





6 12 24 


处 理 时 间 ( 


Treatment time 














图 3 qRT-PCR 检测 miR-2769 靶 基 因 BmETS 不 同 剪接 体 在 20E 处 理 后 不 同时 间 家 蛋 幼虫 脂肪 体 中 的 表达 模式 
Fig. 3 Expression profiles of different splice isoforms of the target gene BmE75 of miR-2769 in larval fat body 





of Bombyx mori at different time after treatment with 20E by qRT-PCR 











A: BmE75B; B; BmE75A; C: BmE75C. 图 中 数据 为 平均 值 上 标准 误 ， 




















] 双 尾 配 对 i 检验 进行 统计 分 析 . Data in the figure are mean + SE, and 


statistics analysis was performed by two tailed, paired t-test. “* P <0.001; “P<0.01; “无 显著 差异 No significant difference (P >0. 05). 


2.3” 双 荧光 素 酶 报告 载体 系统 分 析 miR-2769 与 
BmE75B 的 互 作 

BmE75B 3' UTR 序列 的 PCR 扩 增 结果 如 图 4 
(A) 所 示 , 扩 增产 物 的 分 子 量 大 小 与 理论 一 致 。 扩 
增产 物 连接 psiCheck2 双 效 光 素 酶 报告 基因 载体 后 
的 双 酶 切 鉴定 结果 如 图 4(B) 所 示 , 酶 切 产 物 大 小 与 
理论 分 子 量 相符 。 构 建 好 的 psiCheck2/BmE75B 3' 
UTR 双 葡 光 素 酶 报告 载体 经 测序 验证 正确 。 

将 psiCheck2/ BmE75B 3'UTR 双 荧 光 素 酶 报告 
载体 分 别 与 合成 的 miR-2769 mimics, miR-2769 
mutant 和 miR-2769 NC 共 转 染 HEK293T 细胞 ,以 
psiCheck2/BmE75B 3'UTR 质粒 单 转 为 对 照 , 利 用 双 
荧光 素 酶 报告 基因 检测 试剂 盒 检测 转 染 后 各 组 的 海 
PFET ZR ETE VE at Hk HE Ze ME PE, FAY EL 
(EK HEE RA AT PE. SHR GU S 所 
aR, 5j psiCheck2/BmE75B 3'UTR 质粒 单 转 染 对 照 
相 比 ,只 有 psiCheck2/BmE75B 3' UTR 与 miR-2769 
mimics 共 转 染 时 ,报告 基因 的 荧光 素 酶 相对 活性 显 
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图 4 BmE75B 3'UTR 序列 的 PCR 扩 增 (A) 和 双 荧 光 素 酶 
报告 载体 psiCheck2/BmE75B 3'UTR 的 酶 切 鉴 定 (B) 
Fig. 4 PCR amplification of BmE75B 3'UTR (A) and 

identification (B) of dual luciferase report vector of 
psiCheck2/BmE75B 3'UTR by restriction digestion 

M1; DL2000 Plus DNA Ladder; 1; BmE75B 3'UTR PCR 扩 增 结果 

PCR product of BmE75B 3'UTR; M2: DL2000 DNA Ladder; 2: 

psiCheck2/BmE75B 3'UTR 双 酶 切 鉴定 Identification of psiCheck2/ 

BmE75B 3'UTR with double enzyme digestion. 
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著 降 低 (P <0.05) ,而 psiCheck2/BmE75B 3'UTR 与 
miR-2769 mutant 或 miR-2769 NC 共 转 染 后 的 酶 活 
性 与 对 照 质粒 无 显著 差异 (P >0.05)。 由 此 表明 ， 
miR-2769 能 够 通过 与 靶 基 因 BmE75B 3'UTR 区 的 
互 作 负 调控 BmE75B 在 翻译 水 平 的 表达 。 由 于 
BmE75A 和 BmE75C 的 3'UTR 也 具有 miR-2769 的 
作用 位 点 ,因此 ,推测 miR-2769 也 可 调控 BmE75A 
和 BmE75C 的 表达 。 

2.4 RẸ miR-2769 与 靶 基 因 BmE75 不 同 剪 接 体 




















均一 化 奖 光 素 酶 活性 


Normalized luciferase activity 








的 表达 分 析 1 2 3 4 

利用 qRT-PCR 进一步 明确 miR-2769 AHALA 图 5 psiCheck2/BmE75B 3'UTR 质粒 与 miR-2769 转 染 
BmE75 不 同 剪接 体 在 家 盘 变 态 发 育 不 同 阶段 的 HEK293T 细胞 后 荧光 素 酶 活性 分 析 
表达 是 否 存 在 负 相 关 ,结果 如 图 6 所 示 ,BmE75B 在 Fig. 5 Analysis of luciferase activity in HEK293T cells 
RAIL We HAS E BYE EY) Pe TS tat FA OY 4 ey TT miR- after transfected with psiCheck2/BmE75B 
2769 在 所 有 的 发 育 阶 段 的 表达 量 都 较 低 ; Bm E75A 3'UTR and miR-2769 
和 BmE75C 的 K 达 特 征 也 不 同 于 BmE75B ， BmE75A 1: psiCheck2/BmE75B 3'UTR SA#@Xt HE psiCheck2/BmE75B 3'UTR 
与 miR-2769 相似 在 所 有 发 育 时 期 的 表达 量 都 较 低 vector control; 2: psiCheck2/BmE75B 3'UTR + miR-2769 mimics; 


yi H a 估量 3: psiCheck2/BmE75B 3' UTR + miR-2769 mutant; 4: psiCheck2/ 
BmE75C 则 在 5 龄 末 和 师 早 期 具有 极 高 的 表达 量 ; 马 BmE75B 3'UTR + miR-2769 NC. 图 中 数据 为 平均 值 标准 误 ; 柱 


由 图 中 小 图 可 见 , 虽 然 miR-2769 与 BmE754 的 表达 | 不同 字母 表示 差异 显著 (P<0.05, 单 因素 方差 分 析 , Tukey EE 
量 较 低 ， 但 二 者 的 # 达 仍 存 在 负 相 关 5 由 此 表 明 检验 ) Data in the figure are mean + SE. Different letters above bars 
miR-2769 5 BmE75A , BmE75B 和 BmE75C TER AES indicate significant difference ( P <0.05, one-way ANOVA, Tukey’ s 
































~ 一 ‘LE ` test). 
态 发 育 的 不 同 阶段 均 存在 一 定 程度 的 表达 负 相关 。 ea 

700 ES) e miR-2769 
=- miR-2769 te E æ BmE75A 

6004 = BmE75A a é 

RE 

500 =- BmE75B 5 E 

ni =- BmE75C E 

+) % 400 = 


300 


相对 表 


Relative expression level 


200 


100 


0 
SN AN YN OY DD AN A&W ws 2% 9 96 SM od 9 WS D 
wr Py SPOS SSO oS TEESE S PS CK KY wy 


发 育 阶 段 Developmental stage 














图 6 qRT-PCR 检测 miR-2769 及 其 靶 基 因 BmE75 不 同 剪 接 体 在 家 大 不 同 发 育 阶 段 的 表达 模式 
Fig. 6 Expression profiles of miR-2769 and different splice isoforms of its target gene BmE75 in 














different developmental stages of Bombyx mori by qRT-PCR 
4L2D -5L8D: 分 别 为 4 龄 第 2 天 至 5 龄 第 8 K 2nd day of the 4th instar to the 8th day of the Sth instar, respectively; PP1 — PP9; 分 别 为 师 期 第 1 
-9 天 1st -9th day of the pupal stage, respectively; Al - A3 : 分 别 为 成 虫 第 1 -3 天 1st -3rd day of the adult stage, respectively. 幼虫 去 中 肠 。 图 
中 数据 为 平均 值 上 标准 误 ; 小 图 为 miR-2769 和 BmE75A 调整 坐标 轴 刻 度 后 的 表达 量 。The midgut of larvae was removed. Data in the figure are 


mean + SE. Inset shows the expression of miR-2769 and E75A in adjusted Y axis scale. 


























为 了 进一步 分 析 miR-2769 与 其 靶 基 因 BmE75 ”qRT-PCR 分 别 对 miR-2769 和 BmETS 不 同 剪接 体 进 
不 同 剪 接 体 在 家 看 个 体 中 的 表达 情况 ,首先 以 家 看 ” 行 检测 。 结 果 表明 ,miR-2769 在 5 Wea R FERA 
5 龄 第 3 天 幼虫 不 同 组 织 的 cDNA 为 模板 ,利用 卵巢 中 的 相对 表达 量 较 高 ,其 他 组 织 的 表达 量 相 对 
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BED ; SESE) BmE75SB EAK ON SEA ASP Ae 
TA EAT Be tes, BRK TB R I TR ES A SS EG EF PE 
达 量 无 显著 性 差异 (已 > 0.05) 外 ,其 他 组 织 中 miR- 
2769 与 BmE75B 的 表达 差异 均 为 显著 (P <0.05)。 
而 BmE75A 和 BmE75C 具有 不 同 于 BmE75B 的 表达 
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特征 ,BmE754 在 血 淋 巴 和 卵 偶 中 表达 量 较 高 , 而 
BmE75C 的 整体 表达 量 都 不 高 ,但 miR-2769 与 
BmE75A 和 BmE75C 在 多 个 不 同 组 织 中 的 表达 也 存 
在 显著 差异 (P<0.05)( 图 7)。 
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图 7 qRT-PCR 检测 miR-2769 及 其 靶 基 因 BmE75 不 同 剪 接 体 在 家 看 5 Wey 





不 同 组 织 中 的 表达 模式 


Fig. 7 Expression profiles of miR-2769 and different splice isoforms of its target gene BmE75 in 


1: 头 部 Head; 2; # 


different tissues of the Sth instar larvae of Bombyx mori by qRT-PCR 





皮 Epidermis; 3: 血 淋 





E Hemolymph; 4; 中 肠 Midgut; 5: 马 氏 管 Malpighian tubules; 6: 脂肪 体 Fat body; 7: 丝 腺 Silk gland; 


8: 精 集 Testis; 9: 卵巢 Ovary. 图 中 数据 为 平均 值 + 标准 误 ; 柱 上 不 同 字 母 表示 差异 显著 (P<0.05， 


mean + SE, and different letters above bars indicate significant difference (P <0.05, one-way ANOVA). 














单 因素 方差 分 析 )。Data in the figure are 


3 讨论 


miRNAs 是 一 类 小 分 子 量 的 非 编 码 RNA, 主要 
通过 与 靶 基 因 的 互 作 调控 mRNA 的 降解 和 翻译 抑 
io WREKE miR-2769 在 20E 诱导 早期 的 
旨 肪 体 中 表达 下 调 ( 图 1) , BmE75B 作为 miR-2769 
的 预测 靶 基 因 ,在 20E 诱导 后 家 符 脂 肪 体 中 的 表达 
与 miR-2769 的 表达 相反 。 同 时 20E 4 RA NW E 
中 可 快速 诱导 BmE75A 的 高 表达 , 而 BmE75C 在 
20E 诱导 后 表达 升 高 较 慢 , 且 表 达 水 平 相对 较 低 (图 
3) 。 利 用 NCBI 数据 中 的 Blast 比 对 分 析 表 明 ,家 看 
的 BmE75A, BmE75B 和 BmE75C 与 黑 腹 果 蝇 
DmE75A, DmE75B 和 DmE75C 编码 的 氮 基 酸 序列 
一 致 性 约 为 50% 。 研 究 表明 , 黑 腹 果 蝇 Dm E75A 和 
DmE75B 被 20E 诱导 后 快速 高 表达 , DmE75C 诱导 
后 表达 则 较 慢 且 表 达 量 弱 的 结果 一 致 ( Bialecki et 
al., 2002; Bernardo et al., 2009) ,由 此 可 见 ,E75 不 
同 剪 接 体 在 20E 诱导 后 家 盔 和 果 蝇 中 具有 相似 的 
表达 特征 。 同 时 ,前 人 也 研究 表明 ,20E 也 可 诱导 家 
看 BmE75 基因 的 高 表达 (Matsuoka and Fujiwara, 
2000; Swevers et al., 2002; Shirai et al., 2012) , H. 








BmE75A, BmE75B 和 BmE75C 对 20E 具有 发 育 阶 
段 和 组 织 的 特异 表达 特征 (Li et al., 2015 ) 。 

利用 双 菊 光 素 酶 报告 系统 研究 miR-2769 与 靶 
FAW AVE, FRA PILE miRNAs 调控 丢 基 因 的 
方式 类 似 (Huang et al., 2012; Jiang et al., 2013; 
Ling et al., 2015; Liu et al., 2018, 2019) , miR-2769 
能 够 通过 与 靶 基 因 Bm E75B 3'UTR 区 的 序列 互 作 调 
控 靶 基因 在 翻译 水 平 的 表达 , 由 于 BmE75A, 
BmE75B 和 BmE75C ft} 3' UTR 具有 相同 的 miR- 
2769 互 作 位 点 ,推测 miR-2769 也 可 调控 BmE75A 
和 BmE75C 的 表达 。 结 果 表 明 在 家 不 幼虫 脂肪 体 
中 ,BmE754 和 BmE75B 在 20E 诱导 家 春 脂 肪 体 后 
不 同时 间 的 表达 量 均 与 miR-2769 的 表达 具有 人 负 相 
关 ; 在 正常 家 乔 5 龄 幼虫 的 大 部 分 组 织 中 ,miR-2769 
与 BmE75B 不 同 剪接 体 的 表达 也 都 存在 负 相 关 , 由 
此 表明 ,miR-2769 可 能 通过 与 BmE75 3' UTR 的 互 
作 调 控 其 不 同 剪 接 体 在 不 同 组 织 中 的 表达 。 

不 同 发 育 阶段 和 组 织 中 的 表达 特征 表明 ,miR- 
2769 与 BmE75 的 不 同 剪接 体 在 4 龄 第 2 天 到 成 虫 
不 同 发 育 阶段 和 5 龄 幼虫 不 同 组 织 中 的 表达 均 具 有 
特异 性 表达 特征 。 在 血 淋 巴 中 miR-2769 可 促进 
BmE754 的 表达 ,在 脂肪 体 中 miR-2769 可 促进 
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BmE75B 的 表达 ,而 在 其 他 大 部 分 组 织 中 miR-2769 
抑制 BmETS 不 同 剪接 体 的 表达 (图 7) ,推测 可 能 与 
家 看 不 同 组 织 中 20E 的 滴 度 变化 有 关 。 由 于 家 看 4 
龄 幼虫 到 成 虫 为 整 虫 取样 , 整 虫 样品 中 miR-2769 的 
模板 量 可 能 被 稀释 ,使 miR-2769 FEA EAST] AC IY 
期 的 表达 量 整体 较 低 , 而 BmE75A, BmE75B 和 
BmE75C 在 不 同 阶段 的 表达 量 则 各 不 相同 ,尤其 
BmE75C 在 5 龄 末 和 晴 初 期 表达 量 极 高 ; BmE75B 
在 5 龄 末 、 晴 期 第 1 天 、 肾 期 第 4 天 , 肾 期 第 7 天 和 
羽化 前 的 变态 发 育 时 期 表达 量 较 高 (图 6) ,该 结 
与 家 乔 变 态 发 育 过 程 中 20E 的 滴 度 变化 (Burdette， 
1962 ) 基本 一 致 ; BmE7SA 的 表达 量 则 整体 较 低 ,但 
大 部 分 时 间 点 miR-2769 与 BmETS 不 同 剪接 体 的 表 
达 均 为 负 相 关 。 在 烟草 天 蛾 Manduca sexta 2) HE Wt 
皮 和 化 肾 过 程 中 ,检测 到 E754 和 E75B 的 表达 与 
20E 水 平 的 升 高 一 致 ,但 8754 的 表达 稍 早 于 E75B 
(Zhou et al., 1998), FHF 20E 在 家 在 不同 发 育 阶 
段 及 不 同 组 织 的 滴 度 不 同 ,可 能 使 miR-2769 与 
BmET5 的 不 同 剪接 体 在 不 同 组 织 和 不 同 发 育 时 期 
表现 出 不 同 的 调控 模式 。 另 外 , miR-2769 和 
BmE75B 在 精 集 和 卵 集 中 的 表达 量 非常 高 , 而 
BmE75A H 98 K P HRA miR-2769 相当 (图 
7) ,而 非 负 相关 。 在 家 看 、 烟 草 天 蛾 M. sexta 和 帕 
HE hht Pararge aegeria rP HN rhe Se TK YY miR-2763 
和 miR-989 ,在 卵巢 的 发 育 中 发 挥 作用 (Jagadeeswaran 
et al., 2010; Zhang et al., 2012; Quah et al., 2015) , 
RAE miR-2769 是 否 参与 卵 介 的 发 育 还 有 待 进 一 步 
证 实 。 此 外 ,在 肉 性 果 晶 中 ,DmE754 和 DmE75B 在 
卵 室 发 育 和 凋 亡 过 程 中 发 挥 相 反 的 作用 , DmE754 
可 诱导 哺育 细胞 的 凋 亡 , DmE75B 可 促进 卵 的 发 育 
(Terashima and Bownes, 2006) ,二 者 在 家 乔 卵 发 育 
中 是 否 具 有 相似 的 功能 还 未 见报 道 。 最 近 , Liu 等 
(2018) 发 现 miR-14 也 可 通过 与 靶 基 因 E75 和 ECR- 
B 3'UTR 的 互 作 参 与 家 和 看 幼虫 的 发 育 过 程 ,由 此 推 
WKE BmETS 基因 可 能 还 存在 除 miR-2769 外 的 其 
他 调控 因子 ,参与 其 可 变 剪接 和 功能 调控 。 

综 上 所 述 ,miR-2769 作为 家 看 20E 信号 通路 中 
关键 基因 BmETS 新 的 调控 因子 ,可 通过 与 靶 基 因 
BmE75 3'UTR 区 的 互 作对 BmE75 不 同 剪接 体 进行 
负 调 控 , 同 时 miR-2769 REE BmETS 不 同 剪 接 
体 在 不 同 发 育 时 期 和 不 同 组 织 中 具有 不 同 的 表达 特 
征 。 该 研究 不 仅 丰 富 了 miRNAs 参与 家 看 变态 发 育 
调控 的 功能 研究 ,同时 也 为 家 看 Bm ETS 不 同 剪接 体 
的 调控 机 制 及 其 在 昆虫 变态 发 育 过 程 中 的 功能 研究 














































































































提供 了 新 的 依据 。 
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